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The long-range electromagnetic interaction presents a challenge for numerical computations in 
QCD + QED. In addition to power-law finite volume effects, the standard lattice gauge theory 
approach introduces non-locality through removal of photon zero-momentum modes. The resulting 
finite volume effects must be quantitatively understood; and, to this end, non-relativistic effective 
field theories are an efficient tool, especially in the case of composite particles. Recently an oddity 
related to non-locality of the standard lattice approach was uncovered by the Budapest-Marseille- 
Wuppertal collaboration. Explicit contributions from antiparticles appear to be required so that 
finite volume QED results for a point-like fermion can be reproduced in the effective field theory 
description. We provide transparency for this argument by considering point-like scalars and spinors 
in finite volume QED using the method of regions. For the more germane case of composite parti¬ 
cles, we determine that antiparticle modes contribute to the finite-volume electromagnetic mass of 
composite spinors through terms proportional to the squares of time-like form factors evaluated at 
threshold. We extend existing finite volume calculations to one order higher, which is particularly 
relevant for the electromagnetic mass of light nuclei. Additionally, we verify that the analogous 
finite volume contributions to the nucleon mass in chiral perturbation theory vanish in accordance 
with locality. 

PACS numbers: 12.39.Hg, 13.40.Gp, 13.60.Fz, 14.20.Dh 


I. INTRODUCTION 

Lattice QCD computations have progressed to the level 
of providing a variety of quantitative input for strong in¬ 
teraction phenomenology. Refinements to techniques and 
advances in both algorithms and hardware continue to be 
made, and this has rendered high-precision lattice deter¬ 
mination of some observables in need of QED corrections. 
Such corrections, moreover, are important in address¬ 
ing the light quark mass difference from first principles, 
computing the nucleon mass splitting, and accounting 
for charged hadron scattering at low energies. Each of 
these is a milestone in the further development of lat¬ 
tice QCD. Since the original investigation of QED effects 
on the spectrum of hadrons [T| , there have been consid¬ 
erable improvements in recent lattice calculations 
including an impressive computation of the hadron spec¬ 
trum that dynamically includes all sources of strong and 
electromagnetic isospin breaking Q. 

In tandem with these lattice QCD + QED calcula¬ 
tions, there has also been a focus on addressing theo¬ 
retical aspects of the method and applications beyond 
spectroscopy, for recent progress in these directions, see, 
for example, [lol - [l3l| . In this work, we focus on QED 
corrections to hadron masses at leading order in the hne- 
structure constant, a. Specifically we investigate the per¬ 
nicious finite volume effects that must be quantitatively 
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understood in order to extract physical results. These 
effects arise from altering the long-range behavior of the 
electromagnetic interaction on a torus, and our analy¬ 
sis specializes to the non-local formulation of finite vol¬ 
ume electrodynamics, in which photon zero-momentum 
modes are removedB An ideal framework for these corn- 
stations is that of non-relativistic QED (NRQED) [H- 
Il8l| . which has recently been applied to composite non- 
relativistic particles, such as the nucleon [T^ . 

Within an NRQED framework, finite volume correc¬ 
tions to the electromagnetic masses of hadrons have been 
derived EH- In the case of a point-like spin-half particle, 
results derived in QED disagree with those from NRQED, 
and this is supported by numerical data Quite re¬ 
cently [I^ , it has been argued that the standard NRQED 
description is missing explicit contributions from antipar¬ 
ticles in finite volume on account of the non-locality. In 
the present note, we provide transparency for this argu¬ 
ment by considering the electromagnetic mass of point¬ 
like scalar and point-like spinor particles in QED utilizing 
the method of regions; in the present context, see [lain. 
Additionally we extend these considerations to composite 
scalar and spinor particles, and find there are no mod¬ 
ifications to the finite volume effects from antiparticles 
in the former case, while there are 0{aL~^) effects in 
the latter that depend on squares of time-like form fac¬ 
tors evaluated at the pair-production threshold. Exist- 


^ A local formulation of finite volume QCD + QED using C* 
boundary conditions 0 has quite recently been proposed [T^, 
and shows considerable promise at ameliorating the size of finite 
volume corrections. 
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ing finite volume calculations are trivially extended to 
one order higher. The absence of such corrections is use¬ 
ful to know, for example, in assessing the effect of finite 
volume on the electromagnetic mass of light nuclei. We 
comment briefly on the nucleon-antinucleon contact op¬ 
erators in chiral perturbation theory, however, their finite 
volume contributions to the nucleon mass are shown to 
vanish due to locality of finite volume QCD. 

The organization of this note is as follows. First in 
Sec. m we consider the case of point-like particles in 
QED. For scalars and spinors, we determine the finite- 
volume electromagnetic mass to leading order in a, in¬ 
cluding both power-law and exponential contributions. 
These results are contrasted with those obtained by the 
standard non-relativistic framework, in which the power- 
law contributions should be reproduced. Disagreements 
are linked to antiparticle contributions, as we exhibit by 
employing the (asymptotic) method of regions. Includ¬ 
ing particle-antiparticle contact interactions in the non- 
relativistic effective field theory settles the disagreement. 
These observations are then utilized to extend to the 
more germane case of composite scalars and composite 
spinors in Sec. IIIII Finite volume effects for these com¬ 
posite particles are determined to 0{aL~^). In the scalar 
case, there are no antiparticle contributions at finite vol¬ 
ume, whereas for the spinor case, antiparticle modes show 
up exclusively at 0{aL~^), and are completely depen¬ 
dent on the non-locality introduced by subtracting the 
photon zero mode. By extending the point-like spinor 
case considered in [l^ to that of a composite spinor, 
we find that the required volume effect from anitparti- 
cle modes depends on the square of the time-like Sachs 
magnetic form factor evaluated at threshold. We address 
related concerns in the finite volume computation of the 
nucleon mass using chiral perturbation theory in Sec. lIVI 
Antiparticle modes are argued not to make a contribu¬ 
tion to the volume dependence of the nucleon mass due to 
locality of finite volume QCD. Our findings are summa¬ 
rized in Sec. El details about the one-loop computations 
for point-like particles are contained in Appendix]^ and 
shape coefficients used throughout are evaluated in Ap¬ 
pendix |5] 


II. POINT-LIKE PARTICLES 


To establish the need for antiparticle contributions, 
we evaluate the finite volume corrections to the electro¬ 
magnetic mass of charged, point-like particles. We con¬ 
sider both scalars and spinors in QED to leading order in 
the hne-structure constant. These results are contrasted 
with the non-relativistic approach, and discrepancies in 
the power-law corrections are transparently linked to an¬ 
tiparticle contributions by employing the method of re¬ 
gions. 


_ 

FIG. 1. Electromagnetic mass corrections at 0{a). Circles 
denote the relativistic charge and charge-squared couplings of 
the photon. The latter coupling is only present in the case of 
scalar electrodynamics. 


A. One-Loop Computations 

In all of our computations, we consider the theory 
of QED defined on a three-dimensional spatial torus of 
length L in each direction, with the time direction kept 
infinite. To avoid singularities associated with the long- 
range nature of the electromagnetic interaction, we in¬ 
troduce non-locality by removing the three-momentum 
zero modes of the photon. The effect of this non-locality 
surfaces in finite volume corrections. We work in the par¬ 
ticle’s rest frame, = (M, 0). Results in moving frames 
are not identical because the finite volume breaks Lorentz 
invariance. 


1. Scalar QED 


The relativistic scalar QED Lagrangian density has the 
form 


= ( 1 ) 


where = ((3^ -|- zQeA^)$, along with = 

and e is the magnitude of the electron’s 
charge. To leading order in the fine-structure constant a 
and to all orders in the particle’s Compton wavelength 
M~^, the long-range part of the electromagnetic mass is 
determined by the sunset and tadpole diagrams shown in 
Fig. [TJ Of course there are also short-distance contribu¬ 
tions required which renormalize the particle’s mass. 

Evaluation of the Feynman diagrams is straightfor¬ 
ward. In the rest frame and with the subtraction of the 
photon’s three-momentum zero modes, we arrive at the 
following intermediate expression for the finite volume 
electromagnetic mass shift 


AM$ = —iQ'^e^ 



1 

+ ie 


‘2Kd -f ^0 

2MqQ-\-q^ + ie 2M 


( 2 ) 


Throughout, finite volume corrections to a quantity Ai 
are defined by AAI = M{L) — AI(oo), and for brevity 
we employ the notation 


q n^O 



dq 

(27r)' 


■/(g), 


( 3 ) 


where the hat denotes exclusion of the three-momentum 
zero mode. 
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The finite volume correction shown in Eq. @ evaluates 
to 


provide an intermediate stage of the calculation, at which 
point we have the finite volume mass shift 


AMj, ^ _ 67rQ^a 

M ~ M (27rMT)3/2® ’ ^ 

where, for illustrative purposes, we have retained the 
leading exponential correction, see Appendix for the 
full expression. The universal power-law finite vol¬ 
ume correction appearing above is the same 

for all particles iini, includes terms at 0{aL and 
0{aL~‘^)^ and has the explicit form 

AM(b2) q 2 ^ 

M ~ 2-kML^'^^ 

The parameters Cj above are pure numbers that reflect 
the modes allowed in the spatial geometry; and, for a 
spatial torus, they are defined by 



AM^ = —2iQ^e^ 


dqo 1 M — qQ 

2tt ^ g2 _|_ jg 2Mgo + q^ + ie 

(9) 


Carrying out the computation of the finite volume effect, 
we arrive at 


AMq, AM^^’2) ^ 247rQ^a -ml 

M M (ML)3 ~ (27rMT)3/2® ’ 

( 10 ) 

where the universal finite volume correction is given 
in Eq. ([5]), and only the leading exponential has been 
shown, see Appendix lAl Unlike the scalar case, there is 
a non-vanishing correction beyond the universal power- 
law terms because of the subtraction of photon’s zero- 
momentum mode. The power-law corrections are identi¬ 
cal to those obtained in Q. 


B. Non-Relativistic Particles 


see Appendix IbI for their evaluation. 

In considering non-relativistic effective field theories 
below, only the power-law finite volume dependence can 
be determined. To this end, it is helpful to cast the finite 
volume correction, Eq. in the form 

H- , 

(7) 

where the small dimensionless parameter ^ is defined 
by ^ Higher-order terms neglected above 

all depend exponentially on the volume. Consider¬ 
ing an expansion in the power-law corrections rep¬ 
resent perturbative contributions that can be obtained 
in the effective theory, while exponential corrections are 
non-perturbative and are thus inaccessible in the non- 
relativistic effective theory. 


AAf<j 

M 


= Q^a 




^ 2 exp ( — - 


2. QED 

The interactions of a point-like spin-^ particle 41 in 
QED are described by the Lagrangian density 


Now we wish to address whether the finite volume be¬ 
havior for scalar and spinor particles can similarly be 
obtained using non-relativistic effective theory. To ac¬ 
complish this, we shall transparently derive the non- 
relativistic results using the method of regions [HI ■ 

The advantage of this method is a one-to-one correspon¬ 
dence between regions of integration and the effective de¬ 
grees of freedom. Of course, the disadvantage is that this 
asymptotic method can only be applied in perturbation 
theory, and so will not produce exponential corrections. 
Furthermore, the method does not straightforwardly ex¬ 
tend to composite particles. 

We start from the relativistic expressions for the mass 
shift given in Eqs. m and ®, and inspect the pole struc¬ 
ture of the integrand. In terms of pole structure, both 
expressions are essentially the same. The photon prop¬ 
agator has poles exclusively in the small-energy regime, 
qo ^ |q| ^ M. On the other hand, the propagator of the 
matter field has a pole in this small-energy regime, as well 
as a pole in the regime qo ~ —2M, |q| ^ M, which cor¬ 
responds to the contribution from an intermediate-state 
antiparticle. To exhibit this, we write the propagator of 
the matter field as 


£ = ( 8 ) 

where = (9^ -I- To leading order in the 

fine-structure constant a and to all orders in the parti¬ 
cle’s Compton wavelength the long-range part of 

the electromagnetic mass is determined by the sunset di¬ 
agram shown in Fig. [T] Short-distance contributions that 
renormalize the particle’s mass are also required but not 
depicted. 

As in the scalar case, evaluation of the Feynman dia¬ 
gram in the particle’s rest frame is straightforward. We 


2Mqo + q^ + ie {qo — E++ ie){qo — E-— ie)’ ^ ^ 

with 

E+ = -M+y/APT^=^ + --- , 

E_=-M- = -2M - + • • • , (12) 

where E± correspond to the particle and antiparticle 
states, respectively. 
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In the small-energy region of integration, we can ap¬ 
proximate the antiparticle pole by analytic terms because 
qo <C Inserting the expansion 



and integrating term-by-term clearly leads to an asymp¬ 
totic expansion, but produces the non-relativistic effec¬ 
tive field theory answer. As the poles from the photon 
and particle state are retained, these are the correspond¬ 
ing effective degrees of freedom. For the scalar case, this 
non-relativistic expansion gives rise to the finite volume 
electromagnetic mass shift 



FIG. 2. The fully relativistic sunset diagram reduces to 
the antiparticle tadpole topology of the effective theory in 
the limit of large photon virtuality. The antiparticle tadpoles 
give the contributions in Eqs. (HU and (HI) that are required 
for the effective field theory to reproduce the QED results. 

replacements 

(go - = (po - - 2M - 

(go -E+ + ze)-i = -(2M + E+)-^ + O(po), 

(g2+i£)-i= {4M^-q^)-^+0{po). (16) 


(14) 

where all higher-order power-law corrections, which have 
the form for j = 1,2,-••, are absent due 

to vanishing of the corresponding shape coefficients c_ 2 j, 
see Appendix |B] In this way, we recover the finite volume 
effect that is perturbative in the parameter see Eq. ©. 

For the spinor case, the non-relativistic expansion gives 
us the finite-volume mass shift 

AM^^ STrg^a ^ 

M M 2(ML)3 ^ ^ 

in which the 0{aL~^) contribution is a factor of two 
smaller than the correct result appearin g in Eq. 

This discrepancy has been discussed in [3, [l3| ■ 

Given the discrepancy, we see transparently that the 
finite volume effect from antiparticles cannot be repre¬ 
sented by the series of local terms obtained in Eq. (US. 
Further degrees of freedom are hence required. No¬ 
tice that this analytic antiparticle approximation can be 
found by ignoring issues of convergence, and directly se¬ 
ries expanding the integrands in Eqs. and m as a 
function of M~^ before integration. 


C. Non-Relativistic Antiparticles 

The discrepancy between relativistic and non- 
relativistic results can be resolved by including contri¬ 
butions from antiparticles. We first show that this is 
the case using the method of regions. Then we turn to 
accounting for antiparticle contributions in the effective 
theory. 

Returning to Eqs. m and ®, we investigate the re¬ 
gion of integration where the antiparticle contribution 
is dominant. In this region, we have go ^ —2M and 
I < 7 1 <C M, for which the particle contribution can be ap¬ 
proximated by a tower of analytic terms. The photon 
propagator, moreover, can similarly be approximated by 
analytic terms in this region. Performing a series expan¬ 
sion about po = 9o + 2M r; 0, we have the leading-order 


Notice that -|- 2M vanishes for large M. The phys¬ 
ical interpretation of these approximations is that only 
the antiparticle part of the matter field propagates in 
the loop. The particle part has been replaced by ana¬ 
lytic terms. This is additionally the case for the photon 
propagator, as the photon necessarily has large virtual¬ 
ity, g^ > (2M)^. In the effective field theory description, 
we thus have the external states coupled to virtual an¬ 
tiparticle modes via local operators, see Fig. [2j 

Inserting the above approximations into the scalar elec¬ 
tromagnetic mass shift, Eq. and subsequently ex¬ 
panding the integrand in we arrive at the leading 

contribution from the antiparticle pole 



Evaluation of the integral does not produce a finite 
volume correction at 0{aL~^), and explains why non¬ 
locality has no effect on point-like scalars. 

For the spin-i case, the leading antiparticle contribu¬ 
tion takes the form 





on account of inserting the approximations in Eq. dSD 
into the expression for the mass shift, Eq. ([H|). In turn, 
the missing contribution to the spinor’s finite volume 
electromagnetic mass shift emerges 


AM|"*‘ Stt Q^a ^ 
M ~ ~2{MLf 


(19) 


which is such that AA/^R AM“R = AM^,. 

To capture such missing contributions in an effective 
field theory framework, we thus must retain antiparticle 
degrees of freedom. We begin by considering the scalar 
case, for which we must demonstrate that such contri¬ 
butions are absent. The Lagrangian density including 
the non-relativistic scalar field (^, and the non-relativistic 
anti-scalar field x schematically reads 


C = CptXcj), e) -H Lptix, -e) + Ac((^, x), (20) 
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FIG. 3. Particle-antiparticle annihilation diagrams in the full 
and effective theory. The necessarily large photon virtuality 
reduces to a contact interaction in the effective theory. 


where £c(<^, x) contains contact interactions and is our 
present focus. The non-relativistic fields interact through 
the exchange of photons given the couplings for point¬ 
like particles. These interactions reproduce the photon- 
exchange interaction of the relativistic theory order-by¬ 
order in M~^, with one exception; the annihilation chan¬ 
nel. In this channel, the intermediate-state photon has 
large virtuality which excludes this type of interaction 
from the effective theory. For large virtuality, however, 
the photon propagator contracts to a point, and the an¬ 
nihilation process thus appears in the effective theory in 
the form of contact interactions between particles and an¬ 
tiparticles [i^]. These appear in Cc{(t>,x)- Matching the 
relativistic annihilation diagram to a contact interaction, 
see Fig. O we find 

^c((/>,x) = , (21) 

as the leading contribution in the limit of small momenta. 
In the above expression, we define the Hermitian deriva¬ 
tive iDfj^ by its action AHD^B = A^[iD^B\ — [iD^A^]B. 
Notice from Eq. (HD), that the limit of small momenta 
corresponds to expanding the photon’s virtuality above 
the threshold for pair production. 

Operators of the contact Lagrangian produce contri¬ 
butions to the particle’s self energy through tadpole dia¬ 
grams. The antiparticle tadpoles arise from self contrac¬ 
tion of the antiparticle fields in the contact Lagrangian. 
These are the contributions considered in [l^ for a point¬ 
like spin-half particle. Physically the tadpole diagram 
takes into account the contribution to the relativistic 
sunset diagram in the limit of large photon virtuality, 
see Fig. [2] Computation of the antiparticle tadpole di¬ 
agram in the effective theory generated by the operator 
in Eq. m gives zero on account of the vanishing of c_ 2 . 
Thus the finite volume correction for a point-like scalar 
remains universal in the effective theory. 

There are operators of lower mass dimension in the 
contact Lagrangian, such as 

x) = {4>^x) (xV) , (22) 


that produce finite volume contributions to the electro¬ 
magnetic mass. The contact operator above leads to an 
electromagnetic mass shift of the scalar 


AM| ^ CQ^ 
M 8(ML)3 


( 23 ) 


TABLE I. Low-energy constants for composite hadrons in 
sNRQED and NRQED. Listed are those coefficients required 
in the present study. 


Constant 

sNRQED 

NRQED 

Cf 

- 

Q K 

CD 

<r^ > 

Q + < r| > 

Cm 

> 

< r| > -\k 

cs 

- 

Q + 2/^ 

CAi 



CA2 

fgM^ < > 

— S2'KM^{aE + Pm) 

IQM^ <rl> + 2«:2 
— S27vM^ (oe + Pm) 


which is non-vanishing due to the zero-mode subtraction. 
The coefficient C of this lower-dimensional contact opera¬ 
tor, however, requires matching scalar QED at one-loop, 
for which C ^ a^, and the corresponding finite volume 
correction from antiparticles scales as 0{a^L~^), which 
is beyond our present considerations. 

In the spin-i case, the relevant particle-antiparticle 
contact operator was identified in |12j| . and we summa¬ 
rize their findings. For a non-relativistic spinor field tp 
and its antiparticle Xj the leading contact operator has 
the form 

x) = • (x^c^V’), (24) 

where the coefficient has been determined by matching 
to QED at tree level [l^. The antiparticle tadpole gen¬ 
erated from this spin-dependent contact interaction re¬ 
produces the missing contribution to the electromagnetic 
mass, Eq. m- Having accounted for antiparticle con¬ 
tributions to the finite volume electromagnetic mass of 
point-like particles, we now extend our consideration to 
composite particles using NRQED. 


III. COMPOSITE HADRONS 

Given our demonstration for point-like particles, it is 
natural to extend the discussion to composite hadrons. 
In the context of non-relativistic effective field theories, 
we begin with composite scalars, and then move on to 
composite spinors. In each case, we determine the finite 
volume electromagnetic mass shift to 0{aL~^), which is 
one order beyond that computed in [lll | , and include an¬ 
tiparticle contributions, which are shown to depend on 
squares of time-like electromagnetic form factors evalu¬ 
ated at threshold. 


A. Composite Scalars 

Scalar NRQED (sNRQED) consists of three parts: the 
Lagrangian density for the particle field, that of the 


















antiparticle field, and finally their contact interactions. 
Thus we write the sNRQED Lagrangian density in the 
form 
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Anrqed = C{(l), e) + £(x, -e) + x), (25) 


where the particle and antiparticle fields are described by 
the Lagrangian density 


C{(t),e) = (f)' 


^ [V • E\ 


{D\[VxBY} 2 

-\-IGCm -oTTt-^ ^ 


8M3 
- E'^ 


8M3 


e CA 2 


E^ 


16M3 


+ ecxi 


8M3 

[D^,D-E + E- D] 


{D^iV-E]} , 
+ecx2--h ecx. 


+ie^cx. 


16M4 
{D\{E X By} 
16M4 


16M4 
[V2 V • E] 
16M4 


(26) 


The constants appearing above have been determined 
in [^ . and the ones relevant for the present work are 
given in Table HI 

To the order we work, only the leading contact interac¬ 
tion Lagrangian density is required. To leading order in 
a, the contact operator is given in Eq. m, however, the 
coefficient of this operator must be modified to reflect the 
composite nature of the hadron. Due to the perturbative 
nature of QED, it suffices to consider the one-photon an¬ 
nihilation diagram in Fig. |3] to perform the matching. 
Unlike the point-like case, the annihilation vertex is re¬ 
placed by the full hadronic amplitude, for example 

($(p)<i>t(p')|J^|0) = e{p'-p)^F^{s), (27) 

where s = {p' + p)'^ > 4M^ is the photon’s virtuality, 
and Ej,(s) is the time-like electromagnetic form factor of 
$. Carrying out the computation of the annihilation dia¬ 
gram for a composite scalar and performing the threshold 
expansion leads us to the contact Lagrangian density 

Cc{y, X) = x) • (x^^D , (28) 

where E = F^{s = 4M^) is the form factor at thresh¬ 
old. For a point-like particle, we accordingly recover the 
contact interaction in Eq. (I^Tl) from that in Eq. (05)) . 

With the specification of Eq. (|25|) . we have all ingre¬ 
dients necessary to compute the finite volume electro¬ 
magnetic mass shift of the composite scalar particle to 
0{aL~^). Repeating the calculation at 0{aL~'^), we re¬ 
cover the result derived in [HI- We must augment this 
effective field theory calculation by considering two ad¬ 
ditional classes of diagrams. The first class consists of 
photon-mediated tadpoles, while the second class con¬ 
sists of antiparticle tadpoles. The photon-mediated tad¬ 
poles arise from the particle exchanging a photon with a 
particle or antiparticle loop, see Fig. 01 Both of these con¬ 
tributions are potentially non-vanishing in the effective 


= 0 


FIG. 4. Photon-mediated tadpole diagrams in the effective 
theory. The circles and squares denote the various photon 
couplings present in the NRQED Lagrangian density. For 
each possible square vertex, the sum of particle and antipar¬ 
ticle loops vanishes for a given circle vertex due to charge 
conjugation. 


theory, however, the sum of the two necessarily vanishes 
by charge conjugation invariance. Two photons coupling 
to particle and antiparticle loops are charge conjugation 
even and do not cancel; however, these effects are 0{a‘^) 
and beyond our present considerations. 

The antiparticle tadpole contribution has been deter¬ 
mined above for a point-like particle, and it is trivial to 
modify this result for a composite scalar: 

This contribution is at 0{aL~^), and while there are 
other such terms at this order, they cannot affect the fi¬ 
nite volume mass because the accompanying shape coeffi¬ 
cient vanishes, c 2 = 0. Thus there is no modification to 
the result of [ll| one order beyond what was considered 
in that work. For completeness, we cite the 0{aL~^) 
result here 




-b [ole + Pm) C-i 

+ ^jQ<r^>c., + 0{aL-% (29) 


Notice that in the above expression, there is no con¬ 
tribution from anti-scalar modes. This result is actually 
true to all orders in 0{L~^), and to leading order in a. As 
we have seen, such contributions arise from antiparticle 
tadpole diagrams generated by the contact interactions 
in the effective theory. The coefficients of these operators 
are determined from matching, which requires expanding 
the time-like electromagnetic form factor about thresh¬ 
old, see Eq. (EH). Relative corrections in the threshold 
expansion of the annihilation amplitude always contain 
even powers of the momentum. As a result, the contact 
operators in the effective theory for a composite scalar 
contain an even number of derivatives (which also fol¬ 
lows on account of invariance under parity). Based on 
power counting, the corresponding antiparticle tadpole 
diagrams are thus accompanied by the shape coefficients 
c_ 2 j, for j = 1, 2, • • •. Because all of these shape coeffi¬ 
cients vanish, see Appendix [B] there are no contributions 
from anti-scalar modes in the finite volume electromag¬ 
netic mass of a composite scalar to 0[a). The same 
conclusion does not apply to 0 ( 0 ^) finite volume correc¬ 
tions. 

















7 


B. Composite Spinors 


Now we consider finite volume corrections to the elec¬ 
tromagnetic mass of a composite spin-half particle. By 
retaining the leading contribution from composite spin- 
half antiparticles in NQRED, we derive finite volume cor¬ 
rections valid to 0{aL~^). This requires matching four- 
fermion operators for a composite spinor and generalizes 
the computation of [l^ . 

A composite spinor hadron can be described at low- 
energies by the NRQED Lagrangian density. There are 
three terms necessary: the Lagrangian density for the 
spin-half particle ip, that for the spin-half antiparticle 
X, and the Lagrangian density containing their contact 
interactions. Writing down these three terms, we have 

^NRQED =£^"Hip,e)+C^'‘'>{x,-e) + Ci"Hi^^x), (30) 

where the superscript is used to denote that the La¬ 
grangian density is that relevant for a spinor field. 
For the spin-half particle field, the Lagrangian density 

£(.s) 

can be decomposed into spin-independent and spin- 
dependent terms 

C^^\iP,e) = C{ip,e) + C^'^\ip,e), (31) 


where the spin-independent part C{ip, e) has precisely the 
form given in Eq. (1251) . 

The spin-dependent NRQED Lagrangian density reads 


C^^'^ip, e) = eip^ 


B 


cp- 


2M 


las¬ 


er ■ {D X E - E X D) 


{D^er-B} 

+ CWi- ^7X7^ -^11/2 


-(- Cp'p 


8M3 

{a- D,B- D} 

8^ 


8M2 

DV • BD^ 


ip. 


(32) 


The coefficients of spin-dependent operators that are per¬ 
tinent for our computation are given in Table HI The La¬ 
grangian density for the spin-half antiparticle is related 
by charge conjugation to that of the spin-half particle, 
(X, -e) = £(x, -e) -b (x, -e). 

Unlike the contact interaction in the scalar case, the 
four-fermion contact interaction required in the spinor 


calculation appears at and hence makes a con¬ 

tribution to the electromagnetic mass in finite volume of 
0{L~^). The leading contact operator generated at 0{a) 
is contained in the Lagrangian density 

• (x^'^V')- (33) 

The coefficient of this operator, must be determined 
by matching the hadronic annihilation amplitude with 
that computed in NRQED. As with the scalar case, the 
annihilation process can be computed perturbatively in 
QED using the annihilation diagram shown in Fig. |3l 
The hadron one-photon annihilation vertex must be de¬ 
scribed non-perturbatively by time-like electromagnetic 
form factors, for example 


< '^(.p)'ip(.p')\JtJ.\0) = eu{p) 


Fdshp + ^-^F,is) 


Kp'), 


(34) 

with Qp = p'p Pp, and s = > 4M^. The matching 

computation is an extension of that in [^ . and results 
in the value of the coefficient where 


Gm = Fi{s = 4A/2) +F 2 {s = 4M2), 


(35) 


is the time-like Sachs magnetic form factor evaluated at 
threshold. In the case of a point-like fermion, Gm Q, 
and this ensures the result of [l^ is recovered. 

Having determined the coefficient of the fermion- 
antifermion contact operator, we have all of the necessary 
terms of the NRQED Lagrangian density, Eq. (1551) . to 
perform the finite volume mass calculation to 0{aL~^). 
The result is identical to that derived in m aside from 
the additional term at 0{aL~^) arising from the an¬ 
tifermion tadpole shown in Fig. [21 As with the scalar 
case, the sum of photon-mediated tadpole diagrams 
shown in Fig. |4]vanishes, as do such diagrams with a ki¬ 
netic energy insertion. As a result, there are no 0{aL~'^) 
corrections due to antiparticle modes. One order beyond 
that computed in 0 , we observe that all 0{aL ®) con¬ 
tributions are proportional to the shape coefficient c _2 
which vanishes. For completeness, we give the full result 
to 0{aL~^) accuracy 


AM^ = -b 




47r^ 

u 


(ap + /3m) c_i 


M3L4 


Q ( < r% > -k] - k(Q -b k) 


C-i+0{aL-^). 


(36) 


accompanied by the need to remove loop contributions in 
the effective theory with vanishing antiparticle momen¬ 
tum [H. At higher orders in L however, there are 


Notice that the finite-volume effect from antiparticles in 
the case of spinor particles is due to infrared regular¬ 
ization, namely the subtraction of photon zero modes is 





















no finite volume contributions from antiparticle modes 
at leading order in the fine-structure constant. As in the 
case of a composite scalar, this happens on account of 
parity invariance, which restricts the contact operators 
to have even mass dimension. The corresponding an¬ 
tiparticle tadpoles generated are hence accompanied by 
the shape coefficients C- 2 j, with j = 0, 1, • • •. The only 
non-vanishing antiparticle tadpole on a spatial torus is 
that with j = 0. This contribution has been accounted 
for in Eq. (1551) . 

While the sole 0{a) contribution from antiparticle 
modes is suppressed by two orders relative to the leading 
finite volume correction for charged spinors, the situation 
is different for neutral spinors. The leading finite-volume 
correction for a neutral spinor receives contributions from 
the particle-antiparticle contact interaction. In the case 
of the neutron, for example, we have 

iTrv \ 47r^ 

AM„ = - -JJ + 

(37) 

which is valid to 0{aL~^). As the time-like form fac¬ 
tors of the neutron are not very well known experimen¬ 
tally , careful study of the finite volume effects for the 
neutron electromagnetic mass could, in principle, con¬ 
strain the magnetic form factor at threshold. In practice, 
determining a small finite-volume effect through the elec¬ 
tromagnetic mass of the neutron is quite challenging, as 
the determination of baryon mass splittings already has 
considerable computational demands Q. 

IV. RELATED CHIRAL CONCERNS 

It is instructive to contrast the electromagnetic mass of 
the nucleon with the QCD self-energy of the nucleon. We 
verify that finite volume contributions from the strong 
nucleon-antinucleon annihilation process vanish due to 
the locality of finite volume QCD. 

At low energies, one can use chiral perturbation the¬ 
ory to describe the long-range physics of the nucleon. 
As a result, this effective field theory gives predictions 
for the volume dependence of the nucleon mass in terms 
of the pion Compton wavelength, rn~^. For most lat¬ 
tice QCD applications, one works in the p-regime, where 
itIt^L ^ 1, which ensures that zero modes of the pion 
field do not become strongly coupled [1^. Strong cou¬ 
pling of pion zero modes is quite similar to zero modes in 
finite volume QED, which can be exemplified by giving 
the photon a small mass [1^. The natural question to 
address in light of our work is whether there are annihi¬ 
lation contributions missing from standard chiral pertur¬ 
bation theory analyses of the volume dependence of the 
nucleon mass [ttI . . Nucleon-antinucleon annihilation 
is described by the amplitude depicted in Fig. [5j Un¬ 
like the electromagnetic annihilation process, the strong 
interaction annihilation is not dominated by one-pion ex¬ 
change. The nucleon-pion interactions are weak only at 



FIG. 5. Strong nucleon-antinucleon annihilation. 

low energy, whereas the electromagnetic interaction re¬ 
mains weak well above the pair-production threshold, 
> 4M^. The strong annihilation process can, how¬ 
ever, be described by local operators in chiral perturba¬ 
tion theory. Momentum-independent operators, such as 
the spin/isospin independent contact operator 

^ 

would appear to give the leading contribution from the 
strong annihilation process. In infinite volume, such con¬ 
tact operators lead to an infinite renormalization of the 
nucleon mass, which can be absorbed by local nucleon op¬ 
erators in chiral perturbation theory. At finite volume, 
there is the possibility of a finite counter-term from the 
antinucleon tadpole diagram. Due to momentum inde¬ 
pendence, however, the contributions from all such oper¬ 
ators are proportional to 



The vanishing of this regulated sum depends on the re¬ 
tention of antinucleon zero modes. In the e-regime, where 
niT^L ^ I, a separate analysis is required. 

Antinucleon loops could contribute to the finite volume 
shift of the nucleon mass through higher-dimensional 
contact interactions. Operators of the form 

for example, would lead to finite volume contributions 
oc where the constant of proportionality 

depends the annihilation amplitudes at threshold in the 
possible spin/isospin channels. Such higher-dimensional 
operators with an even number of derivatives, however, 
lead to mode sums involving the shape coefficients C- 2 j, 
which vanish. Thus provided one remains in the p-regime 
of chiral perturbation theory, there are no finite volume 
effects from any antinucleon modes. This is to be ex¬ 
pected from the locality of finite volume QCD. 

V. SUMMARY 

Above we consider the effect of antiparticle contribu¬ 
tions in low-energy effective field theories. Interactions 
of non-relativistic particles and antiparticles can be de¬ 
scribed in NRQED, such as those required in the treat¬ 
ment of positronium. In applying non-relativisitic effec¬ 
tive theories to single hadrons, however, it would seem 
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that such short-distance antiparticle contributions only 
show up at low energies in the form of counter-terms. In 
infinite volume, this is precisely the case as can be demon¬ 
strated from computing annihilation tadpole diagrams. 
Such diagrams, however, can lead to non-vanishing ef¬ 
fects at finite volume. In the non-local formulation of 
finite volume QED, antiparticles make finite contribu¬ 
tions to the electromagnetic mass of hadrons with spin 
at 0{a). 

To make the discussion of these points concrete, we 
consider the case of point-like particles in QED. Antipar¬ 
ticle modes are exhibited in the covariant computation 
using the method of regions. In order for the power-law 
finite volume effect in the non-relativistic description of 
point-like particles to agree with the corresponding QED 
results, we require the addition of particle-antiparticle 
contact interactions. The strength of such interactions 
can be determined from matching conditions, and power 
counting arguments naturally explain why no antiparti¬ 
cle contributions exist in the scalar case, whereas such 
contributions are required in the spinor case. 

We extend the considerations of [l^ to the case of com¬ 
posite particles. In the case of a composite scalar, the 
matching condition is generalized, however, symmetry 
and power counting arguments allow one to conclude that 
anti-scalar modes do not contribute to the finite volume 
effect at leading order in a. Due to the vanishing of the 
shape coefficient c_ 2 , we additionally find that results de¬ 
rived for composite scalars to 0{aL~‘^) in m are valid to 
0{aL~^). For composite spinors, on the other hand, an¬ 
tiparticle contributions modify results at 0{aL~^); and, 
arise from the zero-mode subtraction utilized to regulate 
QED in the infrared. Unlike the point-like spinor case, 
however, the contribution from antiparticle modes is en¬ 
coded in matching the full electromagnetic annihilation 
amplitude. This requires non-perturbative QCD physics, 
and we find that the finite volume effect is sensitive to 
the square of the Sachs magnetic form factor at threshold. 
From precision study of the finite volume dependence of 
the neutron electromagnetic mass, for example, there is 
the amusing possibility to bound the time-like magnetic 
form factor at threshold because such 0[aL~^) effects 
are leading order for neutral particles, see Eq. (l37)) . In 
the course of our investigation, we found that computa¬ 
tions for composite spinors to 0{aL~'^) are valid to one 
order higher due to the identical vanishing of the shape 
coefficient c _2 appearing in 0{aL~^) finite volume ef¬ 
fects. 

Finally we assess the contribution of antiparticle modes 
in finite volume calculations of the nucleon mass. In chi¬ 
ral perturbation theory, the dominant finite volume ef¬ 
fects arise from long-range pions that wind around the 
lattice volume. Nucleons winding around the lattice 
volume are exponentially suppressed relative to these. 
Antinucleon contributions at finite volume are shown to 
be absent due to the retention of pion zero modes, and 
hence the locality of finite volume QCD. 
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Appendix A: Power-Law and Exponential Volume 
Effects 

In this Appendix, we evaluate the one-loop finite vol¬ 
ume corrections to the point-like scalar and spinor elec¬ 
tromagnetic masses. To accomplish this, we follow Q, 
and separate the power-law from exponential behavior in 
ML. Only the leading exponential behavior is quoted in 
the main text. 

In the case of point-like scalars and spinors, we write 
both finite volume corrections in the form 

AM = -f ,5(AM). (Al) 

For point-like scalars, the correction beyond the universal 
one, which is defined as i5(AM$), vanishes exponentially 
fast in the volume. From Eq. m, this correction is given 

by 

x(A M \ f _ 2M -I- go _ 

^ ~ 2M J 2tt 4^ (go + ie){2Mqo + + ie) ’ 

(A2) 

and accordingly does not contain singularities associated 
with propagation of the photon. Furthermore, the zero¬ 
mode contribution to d(AM$) vanishes, and the sum over 
non-zero momentum inodes can be trivially extended to 
all modes. Carrying out the integration over go and utiliz¬ 
ing the Poisson summation formula to perform a Fourier 
transform, we arrive at 

/n2 poo _ 

<5(AM^) = -^ dA{2M-A)J2M^m), 

(A3) 

where Ko(x) is a modified Bessel function of the second 
kind. Here, the absence of = 0 from the winding num¬ 
ber summation ensures that there are no contributions 
which survive the infinite volume limit. The asymp¬ 
totic behavior of the Bessel function, Ko(a:) ^ 
for X 1, leads to the leading exponential correction 
given in Eq. which includes the periodic images with 

For point-like spinors, on the other hand, the finite 
volume correction beyond the universal one, S{AM<j,), 
contains both power-law and exponential dependence on 
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L. The residual power-law dependence is due to non¬ 
locality introduced by the photon zero-mode subtraction. 
From Eq. m, we find the expression 


6{AM^) 


M 


dqo _ M — go _ 

2 tt ^ (go + ie){2Mqo + q'^ + ie) ’ 


(A4) 

which has a form quite similar to the scalar case. Anal¬ 
ogous to that case, we perform the go integration; how¬ 
ever, in order to utilize the Poisson summation formula, 
we must add and subtract the non-vanishing zero mode 
contribution. This is the source of power-law volume de¬ 
pendence in 5(AM^), and the full result takes the form 


the infrared 


s- 


(d) 


_i_l 

o o 


-1 r 


= / ds- 


^r(l) 


•^ 3 ( 0,6 


1 ,(B2) 


with ^ 93 ( 2 , g) is the Jacobi elliptic-theta function, and 
can readily be continued to j < 0 because the quantity 
in brackets vanishes at s = 0 faster than any power of 
s. The 6^"^^ coefficients, on the other hand, potentially 
contain infrared singularities, and are given by 


= lim 


7r-5r(§) 



(B3) 


£7/0 

(A5) 

The asymptotic expansion of the modified Bessel function 
yields the leading exponential correction to the point-like 
spinor mass presented in Eq. (nnj. 


where the last term is the subtraction of the photon zero 
mode, and the incomplete gamma-function Bj (/r) has the 
definition 


/ OO 


'1 

r 




00 

E 

n—O 


1 

n\ 


-V)' 


n-k I 


(B4) 


Appendix B: Shape Coefficients with Zero-Mode 
Subtraction 


From explicit evaluation, we see the infrared singular be¬ 
havior of the coefficients 


In this Appendix, we evaluate the shape coefficients, 
Cj, which are formally defined by 



These mode number sums arise in finite volume compu¬ 
tations of QED with photon zero-mode subtraction, and 
have been employed in the main text. The sum and inte¬ 
gral are each divergent in the ultraviolet, however, their 
difference is finite in a large class of ultraviolet regulariza¬ 
tion schemes. For this reason, we can treat the ultraviolet 
regularization as implicit. To evaluate these coefficients, 
we add a small mass term to regulate the behavior in 
the infrared, and work in d spatial dimensions. Thus we 
consider 


(d) 

S- 


lim 
(1 —^0 



f d^n 

^0 

[n2-k^2p/2 


(Bl) 


with Cj = Cj the desired quantities of interest. Fol¬ 
lowing we consider the coefficients as functions of 

j > 0, and employ analytic continuation for j < 0. 

We write these coefficients in terms of two contribu¬ 
tions -k . The are all non-singular in 


^ The same conclusions are reached using Liischer’s zeta-function, 
as rigorously detailed in [ 30 II . 




j > d 
j = d 


(B5) 


For any charged particle, the largest possible value for 
j entering the effective theory computation of the self 
energy is j = 2, and arises at leading order in the non- 
relativistic expansion. Not surprisingly, we find that the 
photon zero-mode subtraction eliminates all infrared sin¬ 
gularities provided that we work in d > 2 spatial dimen¬ 
sions. 

For all values of j such that j < d, the coefficients 
are hnite [ 2 ^ . Despite the many possible special 

cases depending on the sign of j, and the evenness or 
oddness of both j and d, we find that the coefficients can 
be compactly written in the form 


Ad) _ 2d 

r(|)j(d-i)’ 


(B6) 


with a limit implied for j = 0, — 2, — 4, • • •. In particular, 
the special case j = 0 has the corresponding value b^f^ = 
— 1 in all spatial dimensions d. Combining with the non¬ 
singular piece, = 0, we have 


= -1. (B7) 

Finally both and bj‘^^ contributions vanish at nega¬ 
tive even integers. This leads us to 

cf =0, for j = -2,-4,..., (B8) 






















11 


in all spatial dimensions d. In three spatial dimen¬ 
sions, one has a special relation between shape coeffi¬ 
cients, In the main text, we have made 

use of the additional values: c\ = —2.83729, and 
= -0.266596. 

As many of the conclusions drawn in the main text 
depend on the values derived in Eq. (IB8I) . we show the 
vanishing of these coefficients alternatively by using an 
explicit Gaussian ultraviolet regulator. Introducing the 
Gaussian regularization and without a mass term, we 
have 



Notice that the mode summation has been extended to 
include the zero mode, because the summand vanishes 
there. Making use of the smoothness of the Gaussian we 
can write 


c(t=limf-^ 


-2i 


A^-O 


d\ 




— \rnr 





0, (BIO) 


where we have used limA-j-o A“e = 0, for any a and 


for /3 > 0. The value of can also be obtained from 
Gaussian regularization in the limit that j —>■ 0. Notice 
the above expression produces the quantity Cq'^^ -|- 1 = 0 
because the zero mode has been included in the summa¬ 
tion. 

Separation of the ultraviolet behavior and volume de¬ 
pendence holds for a large class of regulators, such as 
dimensional regularization and smooth cutoffs. Sharp 
ultraviolet cutoffs, by contrast, tend to produce power 
divergences in finite volume corrections which spoil the 
separation of ultraviolet from infrared physics [29|. In 
hard-cutoff regularization, for example, we have the con¬ 
flicting result 


s(i) 


lim 

A^oo 


A 

- 1 + E 1 

n——A 



0 ^- 1 , 


(Bll) 

where the first term is the explicit subtraction of the 
zero mode. The hard cutoff can be ruled out on physi¬ 
cal grounds, however, as further shape coefficient exhibit 
ultraviolet singularities, e.g.. 


c^o = lim 
A—^■OO 

(B12) 

This behavior is physically unacceptable because the in¬ 
finite volume limit of quantities can no longer be taken 
independent of the ultraviolet regularization. 


.n——A 


dn lA 


= lim = oo. 
A—^oo 
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